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Abstract 

Performance  equations  that  describe  the  dependence  of  cell  potential  on  current  density  for  polymer  electrolyte  fuel  cells  (PEFCs) 
have  been  developed  in  algebraic  form.  The  equations  are  derived  from  the  reduction  of  a  one-dimensional  model  that  takes  into  account, 
in  detail,  the  limitations  of  reactant  transport,  proton  migration  and  electron  conduction,  as  well  as  electrochemical  reactions  within  the 
cathode,  anode  and  membrane  electrolyte.  Reduction  of  the  one-dimensional  model  is  implemented  by  approximating  the  profiles  of 
reactant  concentration  and  ionomer  potential  with  appropriate  functions  and  by  lumping  the  overall  reaction  rates  at  the  reaction  centres 
of  the  catalyst  layers.  Since  the  performance  equations  originate  from  a  mechanistic  one-dimensional  model,  all  parameters  appearing  in 
the  equations  have  a  precise  physical  significance.  In  addition,  individual  potential  losses  caused  by  the  various  limiting  processes  can  be 
clearly  quantified  in  the  equation.  Particularly,  the  equations  for  potential  losses  relevant  to  the  anode  limiting  processes  are  first  revealed 
by  the  present  work.  Thus,  they  can  be  used  as  a  diagnostic  tool  for  PEFC  performance.  Computational  results  show  that  the  performance 
equations  agree  well  with  the  original  one-dimensional  model  over  an  extensive  parameter  range.  The  present  performance  equations  allow 
for  an  efficient  evaluation  of  PEFC  performance  since  the  complexities  of  the  one-dimensional  model  and  the  procedures  for  the  numerical 
solutions  are  completely  avoided.  As  compared  with  previously  developed  performance  equations  dealing  with  PEFC  cathodes  only,  the 
present  equations  are  able  to  provide  accurate  interpretation  on  the  polarization  behaviour  of  a  complete  PEFC. 
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1.  Introduction 

Fuel  cells  are  power-generation  devices  that  directly  con¬ 
vert  the  chemical  energy  of  fuels  and  oxygen  into  electricity 
through  electrochemical  reactions.  Among  all  the  various 
types  of  fuel  cells,  the  polymer  electrolyte  fuel  cell  (PEFC) 
is  the  one  characterized  by  using  a  proton-conductive  mem¬ 
brane  as  its  electrolyte  layer.  Due  to  the  advantages  of  low 
emissions,  low  operating  temperature,  and  short  transient 
period  to  steady-state  operation,  PEFCs  are  promising  can¬ 
didates  for  primary  power  devices  for  electric  vehicles.  Up 
to  the  present  stage,  however,  the  commercialization  of  PE¬ 
FCs  is  still  restricted  by  prohibitively  high  cost.  Reduction 
of  the  cost  of  PEFCs  is  achievable  through  the  develop¬ 
ment  of  improved  components  for  better  performance,  such 
as  catalysts  of  high  activity  for  the  oxygen  reduction  reac¬ 
tion,  membranes  of  low  resistance  for  proton  migration  and 
gas  diffusers  with  minimal  water  accumulation.  On  the  other 
hand,  mathematical  modelling  provides  deeper  insights  into 
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the  transport  and  electrochemical  processes  taking  place  in 
PEFCs,  and  thus  can  be  used  as  an  efficient  means  for  deter¬ 
mining  the  most  appropriate  operating  conditions  and  opti¬ 
mal  design  parameters. 

Mechanistic  models  of  PEFCs  are  formulated  based  on 
phenomenological  transport,  The  kinetics  of  electrochemical 
reactions,  and  mass  and  energy  conservation  laws.  Detailed 
consideration  of  the  geometric  configuration  of  PEFCs  leads 
to  two-dimensional  and  three-dimensional  models  [1-10]. 
In  PEFCs,  the  cell  thickness  is  orders  of  magnitude  smaller 
than  the  other  dimensions.  Thus,  although  one-dimensional 
models  are  quite  simple,  they  are  still  able  to  provide  a  good 
account  of  polarization  characteristics  under  various  oper¬ 
ating  current  densities.  A  wide  variety  of  one-dimensional 
models  have  been  proposed  [11-18]. 

The  cell  potential  of  a  PEFC  can  be  evaluated  as  the  equi¬ 
librium  value  subtracted  by  the  potential  losses  caused  by 
the  limitations  of  reactant  transport,  proton  migration,  elec¬ 
tron  conduction  and  electrochemical  kinetics.  Based  on  this 
observation,  empirical  models  are  formulated  as  a  summa¬ 
tion  of  a  few  terms  that  provide  quantitative  estimations 
of  the  equilibrium  potential  and  the  overpotentials  perti¬ 
nent  to  the  related  limiting  processes.  One  of  the  simplest 
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Nomenclature 

Notation 

effective  platinum  surface  area  per  unit 
volume  (cm2  cm-3) 

b 

dimensionless  coordinate  of  anode- 
diffuser  face 

ch2 

hydrogen  concentration  (mole cm-3) 

dc 

gas-diffuser  thickness  (pan) 

effective  binary  diffusion  coefficient  for 
i  and  j  species  (cm2  s-1) 

< 

effective  diffusivity  of  dissolved 
hydrogen  in  the  catalyst  layer  (cm2  s-1) 

Deff 

uo2 

effective  diffusivity  of  dissolved  oxygen 
in  the  catalyst  layer  (cm2  s-1) 

/a?  fc 

parameter  (V-1),/a  defined  by  Eq.  (7c), 
fc=4acF/RT 

F 

Faraday’s  constant  (96487  C  per 
equivalent) 

Hr2 

Henry’s  constant  for  hydrogen 
(atm.  cm3  per  mole) 

Ho2 

Henry’s  constant  for  oxygen 
(atm.  cm3  per  mole) 

i°  i° 

exchange  current  density  at  the 
reference  condition  (A  cm-2) 

/ 

current  density  (A  cm-2) 

I a » 

characteristic  current  density  (A  cm-2), 

7a  defined  by  Eq.  (If), 

h  =  4FPcDfl_/HoA 

him,  a 

anode  limiting  current  density  (A  cm-2) 

neff  neff 
Ka  >  * c 

effective  protonic  conductivity  for 
ionomer  phase  (S  cm-1) 

protonic  conductivity  for  membrane 
(S  cm-1) 

Ni 

mole  flux  of  species  i  (mole  cm-2  s-1) 

Pa,  Pc 

total  pressure  (atm.) 

R 

the  universal  gas  constant 
(8.314  J  per  mol  K) 

n,  r2 

ratio  of  diffusivities,  defined  by 

Eqs.  (19a)  and  (19b) 

T 

temperature  (K) 

V 

catalyst  potential  (V) 

w 

equilibrium  potential  of  single  cell  (V) 

Va,  Vc 

potential  of  anode  and  cathode  (V) 

v°  v° 

y  a  ’  v  c 

equilibrium  potential  of  anode  and 
cathode  (V) 

W 

dimensionless  position  where  hydrogen 
depletion  occurs 

Xi 

mole  fraction  of  species  i 

A 

Xi  in  the  bulk  flow 

xh2 

vh2  at  the  anode  catalyst 
layer/membrane  interface 

XU2 

vh2  at  the  anode  catalyst 
layer/anode-diffuser  interface 

xo2 

xo2  at  the  cathode  catalyst 
layer/cathode-diffuser  interface 

z 

coordinate  perpendicular  to  the  face  of  the 
anode  diffuser  (|xm) 

Greek  letters 

Oic 

electrode  transfer  coefficient 

Pa,  Pc 

parameter  (S  cm-2),  pa  defined  by 

Eq.  (7b),  pc  =  kf/Sc 

Pm 

parameter  (S  cm-2),  defined  as  km/8m 

8  a ,  8C 

catalyst  layer  thickness  (p,m) 

effective  porosity  of  gas  diffuser 

0 

ionomer  potential  (V) 

^ 'Pa.  9  (Pc 

parameter,  <pa  defined  by  Eq.  (7a), 
cpc  =  Aci"Ho2S2c/4FPcDf 

CTeff  CTeff 
ua  ’  c 

effective  electronic  conductivity  for  gas 
diffuser  (S  cm-1) 

dimensionless  coordinate,  defined  by 

Eq.  (7e) 

Subscripts 

a 

anode 

c 

cathode 

m 

membrane 

is  the  model  proposed  by  Srinivasan  et  al.  [19],  who  sug¬ 
gested  an  equation  for  the  conditions  of  activation  and  ohmic 
control.  Such  a  simple  equation  is,  however,  not  only  un¬ 
able  to  interpret  experimental  data  at  high  current  densi¬ 
ties  but  is  also  incapable  of  elucidating  the  effects  of  most 
operating  parameters.  This  equation  was  further  modified 
by  Kim  et  al.  [20]  through  adding  an  exponential  term  to 
compensate  for  the  potential  loss  at  high  current  densities. 
Squadrito  et  al.  [21]  proposed  an  empirical  equation,  which 
includes  a  logarithmic  function  to  account  for  the  effects  of 
mass-transport  limitations.  Since  all  these  empirical  equa¬ 
tions  are  abundant  with  fitting  coefficients,  they  can  only 
be  used  as  a  data-fitting  tool.  Amphlett  and  his  coworkers 
[22-24]  derived  the  performance  equations  of  PEFCs  in  a 
semi-empirical  way,  in  which  both  mechanistic  and  fitting 
coefficients  are  present  in  the  equations.  Pisani  et  al.  [25]  pre¬ 
sented  a  semi-empirical  derivation  of  performance  equations 
with  the  goal  of  having  the  largest  number  of  mechanistic 
coefficients.  Starting  from  a  macro-homogeneous  model  for 
the  cathode  catalyst  layer  of  PEFCs,  Eikerling  and  Korny- 
shev  [26]  made  efforts  to  formulate  analytical  expressions 
for  the  solutions  of  the  model  equations  in  order  to  ratio¬ 
nalize  the  explicit  effects  of  model  parameters  on  cathode 
performance.  Nevertheless,  only  solutions  for  a  few  limiting 
conditions  were  obtained. 

In  a  previous  work  [27],  a  mechanistic  approach  was 
proposed  to  formulate  the  performance  equations  of  PEFC 
cathodes.  In  this  approach,  the  performance  equations  are 
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reduced  from  a  one-dimensional  model  of  PEFC  cathodes 
by  assuming  appropriate  profiles  for  oxygen  concentration, 
ionomer  and  catalyst  potentials.  The  derived  performance 
equations  are  characterized  in  an  explicit  algebraic  form  as 
empirical  models  while  keeping  the  same  level  of  accuracy 
of  the  original  one-dimensional  model.  In  addition,  all  pa¬ 
rameters  appearing  in  the  equations  are  endowed  with  a  pre¬ 
cise  physical  significance.  Moreover,  all  contributive  terms 
of  potential  losses  caused  by  limiting  transport  processes  and 
electrochemical  kinetics  can  be  quantitatively  estimated.  The 
previous  performance  equations  only  deal  with  the  potential 
losses  pertinent  to  the  cathode  of  PEFCs.  At  medium  and 
high  current  densities,  the  overpotentials  associated  with  the 
anode  and  membrane  electrolyte  might  also  contribute  sig¬ 
nificantly  to  the  overall  potential  loss  of  PEFCs.  Thus,  it  is 
also  desirable  to  develop  performance  equations  that  include 
the  potential  losses  of  the  cathode,  anode  and  membrane 
electrolyte  so  as  to  elucidate  the  polarization  behaviour  of 
a  complete  PEFC.  In  order  to  achieve  this  goal,  attention  is 
first  directed  to  the  transport  and  electrochemical  processes 
occurring  within  the  PEFC  anodes  and  equations  for  anode 
performance  are  formulated.  The  performance  equations  for 
a  complete  PEFC  are  then  derived  by  combining  the  poten¬ 
tial  losses  of  its  cathode,  anode  and  membrane  electrolyte. 
Because  the  accuracy  of  the  cathode  performance  equations 
has  been  explored  previously,  it  is  only  necessary  to  inves¬ 
tigate  the  approximation  errors  of  the  anode.  In  addition, 
the  characteristics  of  the  performance  equations  of  a  single 
cell  are  discussed  and  the  potential  losses  stemming  from 
its  cathode,  anode  and  membrane  electrolyte  are  analyzed. 


2.  One-dimensional  model 

The  one-dimensional  model  considered  in  the  present 
work  is  a  steady- state  and  isothermal  one.  The  present  work 
extends  the  previous  results  for  cathode  performance  to  a 
single  cell,  which  consists  of  a  cathode,  an  anode  and  a 
membrane  electrolyte.  Except  that  the  potential  loss  caused 
by  electron  conduction  in  the  cathode  catalyst  layer  is  ne¬ 
glected,  the  assumptions  employed  previously  for  formulat¬ 
ing  the  cathode  model  are  also  applied  and  will  not  be  re¬ 
peated  here.  Only  the  postulations  concerning  the  anode  and 
membrane  electrolyte  are  listed  below. 

(i)  The  Stefan-Maxwell  equations  are  employed  for 
multi-component  gas  transport  in  the  anode  diffuser. 
It  is  assumed  that  a  hydrogen-rich  stream  coming 
from  methanol  steam-reforming  is  used  as  an  anode 
feed,  which  contains  hydrogen,  carbon  dioxide  and 
saturated  water  vapor.  The  effective  binary  gas  diffu- 
sivity  is  described  by  the  Bruggeman  expression  [12], 
which  includes  the  effects  of  porosity  and  tortuosity. 

(ii)  Water  vapour  remains  saturated  within  the  gas  diffuser 
and  the  mole  flux  of  carbon  dioxide  is  zero  due  to  its 
inertness. 


(iii)  Catalyst  particles  and  conductive  ionomers  are  homo¬ 
geneously  mixed  in  the  anode  catalyst  layer  and  the 
macro-homogeneous  assumption  is  applied;  it  implies 
constant  physical  and  chemical  properties  within  the 
catalyst  layer. 

(iv)  The  hydrogen  oxidation  reaction  follows  the  mecha¬ 
nism 

H2  +  2S^2H-S  (1) 

H-S  e"  +  H+  +  S  (2) 

where  S  denotes  an  active  site  on  the  catalyst  sur¬ 
face.  The  above  mechanism  yields  a  half-order  de¬ 
pendence  of  reaction  rate  on  hydrogen  concentration. 
The  Butler- Volmer  equation  [28,29]  is  applied  for  the 
rate  expression. 

(v)  Henry’s  law  holds  for  the  phase  equilibrium  of  hydro¬ 
gen  at  the  interface  between  the  anode  diffuser  and 
the  anode  catalyst  layer. 

(vi)  The  catalyst  layer  is  nearly  flooded  and  Fick’ s  law  is 
applied  for  the  transport  of  dissolved  hydrogen  within 
it.  The  membrane  is  impermeable  to  dissolved  hydro¬ 
gen. 

(vii)  Potential  loss  due  to  electron  conduction  in  the  anode 
catalyst  layer  is  considered  to  be  negligible. 

(viii)  The  membrane  is  fully  hydrated  and  its  ionic  conduc¬ 
tivity  is  a  constant. 

Based  on  the  assumptions  stated  above,  the  equations 
based  on  mass  conservation,  electro-neutrality  and  Ohm’s 
law  can  be  expressed  in  partially  dimensionless  form.  In  the 
anode  gas  diffuser  (1  <  f  <  b): 


P&  dvH2 

RT8a  d£ 

=  “(I  -  *w,a  -  *H2)Afit2 


X 


1 


_DH2-C02  Xh2^w-C02  +  (1  Xw,a  Xh2)jDH2-w_ 


+ 


w,a 


iCff 


eff 


(3) 


In  the  anode  catalyst  layer  (0  <  £  <  1): 


d2XH2 

d£2 


+  'Pa{exp[/a(Va°  ~  V  +  </>)] 


—  exp  [-MV°-V  +  mx% 
d2XH2  fiad  2<p 


1/2 


=  0 


dC 


+ 


h 


=  0 


(5) 

(6) 


The  model  parameters  appearing  in  the  above  expressions 
are  defined  by: 


fed  - 


2  FP&Df2 


(7a) 
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kf 

Pa  ~  „ 

^a 

(7b) 

Oia  F 

fa  = 

J  RT 

(7c) 

II 

1 

(7d) 

//r2Ch2 

*H2  = 

r  a 

(7e) 

j  _  2^a< 

Hu28a 

(7f) 

where  DfJ_-  is  an  effective  binary  diffusivity  in  the  porous 
medium  for  i  and  j  species;  A;  the  mole  flux  of  species  i\ 
Xi  the  mole  fraction;  a a  the  anode  transfer  coefficient;  Aa 
the  effective  platinum  surface  area  per  unit  volume;  i®  the 
anode  exchange  current  density  at  the  reference  condition; 
T  the  temperature;  R  the  universal  gas  constant;  Pa  the  an¬ 
ode  pressure;  ch2  the  dissolved  hydrogen  concentration; 
the  anode  equilibrium  potential;  V  is  the  catalyst  potential; 
0  is  the  ionomer  potential;  Sa  is  the  thickness  of  the  an¬ 
ode  catalyst-layer;  Hn2  is  the  Henry’s  constant  for  gaseous 
hydrogen  and  its  dissolved  form  in  the  ionomer  phase.  In 
the  above  expressions,  vh2  stands  for  the  hydrogen  mole 
fraction  in  the  gas  phase  of  the  diffuser  but  represents  the 
dimensionless  concentration  of  dissolved  hydrogen  in  the 
ionomer  phase  within  the  anode  catalyst  layer,  as  defined 
by  Eq.  (7e).  The  same  symbol  is  used  in  both  regions  be¬ 
cause,  by  such  definitions,  the  xu2  profiles  are  continuous 
across  their  boundary.  It  is  also  assumed  that  the  transfer 
coefficients  for  anodic  and  cathodic  reactions  at  the  anode 
are  equal;  they  are  denoted  as  aa  in  the  above  equations. 

At  the  anode  diffuser  face  (f  =  b ),  Eq.  (4)  is  applied. 
Besides,  the  hydrogen  mole  fraction  is  equal  to  its  value  in 
the  bulk  flow.  Therefore, 

xh2  =  *H2  (8) 

At  the  interface  between  the  anode  diffuser  and  the  catalyst 
layer  (f  =  1),  in  addition  to  Eq.  (2),  other  conditions  are 
written  as 


vh2  (catalyst  layer)  =  xu2  (gas  diffuser) 

(9) 

Ia  dVHo 

2 F  d£  (catalyst  layer)  =  Nu2 

(10) 

d0 

-  0 

(11) 

dC 

The  catalyst  potential  is  related  to  the  anode  potential,  de¬ 
noted  as  Va,  by 

2FNW 

V  (catalyst  layer)  =  Va  H -  (12) 

Aa 

and 

aeff 

Aa  =  -f-  (13) 

da 


where  crff  is  the  effective  electric  conductivity  of  the  anode 
diffuser  and  da  is  its  thickness.  At  the  interface  between 
the  membrane  and  the  anode  catalyst-layer  (f  =  0),  the 
boundary  conditions  are  formulated  as 


0  =  0 


(14) 

(15) 


The  cathode  equations  are  similar  to  those  of  the  anode. 
Details  have  been  presented  elsewhere  [27]. 

The  method  of  collocation  on  finite  elements  based  on 
cubic  B- spline  interpolation  was  employed  for  the  solutions 
of  the  model  equations  [30].  The  equations  of  anode  and 
cathode  were  solved  separately.  First,  the  cathode  equations 
were  iterated  with  a  fixed  cathode  potential.  After  a  solution 
was  obtained,  the  current  density  was  evaluated,  which  was 
then  employed  and  fixed  during  the  iterations  of  the  anode 
equations.  The  details  of  numerical  procedures  can  be  found 
in  the  previous  work  [27].  The  output  potential  of  a  single 
cell  was  calculated  as  its  equilibrium  value  deducted  by 
the  potential  losses  of  its  anode,  cathode  and  membrane 
electrolyte. 


3.  Anode  performance  equations 


It  is  assumed  that  the  dimensionless  hydrogen  concentra¬ 
tion  profile  within  the  anode  catalyst  layer  be  approximated 
by: 

*h2  =  {[(42)1/2  -  (42)1/2k2  +  (42)1/2)2  (i6) 

Where  xsKi  represents  the  value  of  xu2  at  the  interface  be¬ 
tween  the  anode  catalyst-layer  and  the  anode  diffuser,  and 
xlU2  is  the  value  at  the  interface  of  the  anode  catalyst-layer 
and  the  membrane.  Eq.  (16)  automatically  satisfies  the 
boundary  conditions  required  by  the  vh2  profile  at  f  =  0 
and  1.  Because  the  current  density  is  equivalent  to  the  flux 
of  dissolved  hydrogen  at  the  interface  between  the  anode 
catalyst-layer  and  the  anode  diffuser  multiplied  by  —2 F, 
xlK2  can  be  expressed  as: 

2 

(42)_1  (17> 

where  I  stands  for  the  current  density.  Eq.  (1)  can  be  solved 
analytically,  thus  one  has 


xh2  - 


xh2  ~ 


44 


I  =  -2  FNH?  =  - 


2FraDf_CXh 


RTd 


a 


x  In 


nxu  +  +  xwa 


*w,a/(*w,a+ri(l  -*w,a)) 


X 


T\X  h2  +  r2XSQQ2  +  -Tw,a 
b  \  I'l  (1— •Kw,a)/(-Kw,a+l'l  (1— -Kw,a)) 

co2 


co2 


(18) 
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where: 


n  = 


r2  = 


Deff 

UC()2-w 

Deff 

^h2-co2 

Deff 

^H2-w 

neff 

^h2-co2 


(19a) 


(19b) 


The  anode  limiting  current  (/iim,a)  can  be  obtained  by  setting 
v pj  =  0  and  xsCq2  =  1  —  Jcw,a  in  Eq.  (18).  Similar  to  the 
cathode  [27],  an  explicit  approximation  is  used  to  calculate 
x jj  instead  of  Eq.  (18),  which  has  the  form: 


X  _  vb  _  b 

xH2  —  xh2  xco2 


im,a 


-  1 


(20) 


Eq.  (20)  provides  a  convenient  means  for  evaluating  42 
since  it  skips  the  iterative  procedures  required  by  the  original 
analytical  solution.  It  yields  exact  values  at  I  =  0  and  /iim,a 
and  also  quite  good  accuracy  for  other  intermediate  values. 

The  rate  of  the  hydrogen  oxidation  reaction  is  lumped  at 
a  reaction  centre,  which  is  defined  by 

Jo  24x^4  -  37 


^  x  _ 


JoVW1?  48*H27a  87 


(21) 


Based  on  this  approximation,  it  is  further  postulated  that 
the  ionomer  potential  is  flat  from  the  interface  between  the 
anode  diffuser  and  the  anode  catalyst-layer  to  the  reaction 
centre  due  to  the  absence  of  proton  migration.  It  decreases 
linearly  from  the  reaction  centre  to  the  interface  between 
the  anode  catalyst-layer  and  the  membrane  due  to  a  constant 
protonic  current.  Therefore, 


0  = 
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for  f  f  <  1 

(22) 

According  to  the  above  approximations,  the  current  density 
can  be  evaluated  as: 


I  =  <Pa4{exp[-/a(Va0  -  V  +  (/>*)] 
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where  0*  denotes  the  ionomer  potential  at  the  reaction  center 
and  jcy  2  is  defined  as 
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By  inserting  Eqs.  (12),  (22)  and  (24)  into  Eq.  (23),  the  anode 
potential  can  be  solved  and  expressed  as 


*4  =  Va°  + 


24x^4-37 
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The  overall  potential  loss  in  a  PEFC  anode  is  a  linear  com¬ 
bination  of  individual  contributive  terms  stemming  from  the 
limitations  of  the  involved  transport  processes  and  the  elec¬ 
trochemical  reaction,  which  include  diffusion  overpotential 
and  ohmic  loss  for  electron  conduction  within  the  anode 
diffuser,  diffusion  overpotential  and  ohmic  loss  for  proton 
migration  within  the  anode  catalyst  layer,  and  activation 
overpotential  for  the  hydrogen  oxidation  reaction.  Recall 
that  the  ohmic  potential  loss  caused  by  electron  conduction 
within  the  anode  catalyst  layer  is  not  considered  here  due 
to  its  negligibly  small  value.  The  activation  overpotential 
is  referred  to  as  the  potential  loss  required  to  generate  the 
required  current  density  under  the  condition  without  any 
transport  limitations  and  ohmic  losses.  Thus,  the  activation 
overpotential  for  the  hydrogen  oxidation  reaction  is 


(26) 


The  ohmic  loss  caused  by  proton  migration  within  the 
anode  catalyst  layer  is  estimated  as  the  difference  of  ionomer 
potential  at  the  interface  between  the  membrane  and  the 
anode  catalyst-layer  and  at  the  interface  between  the  anode 
catalyst-layer  and  the  anode  diffuser,  which  has  the  form: 


ohmic  loss  (protonic)  of  anode  catalyst  layer 
2442/a  -  37 
48427a  -  87 

The  diffusion  overpotential  of  the  anode  diffuser  arises 
from  the  decrease  in  hydrogen  concentration  from  the  dif¬ 
fuser  face  to  the  interface  between  the  anode  catalyst-layer 
and  the  anode  diffuser,  which  is  expressed  as 


diffusion  overpotential  of  anode  diffuser 

1  xh2(4^4  +  y]l2x  h2  +  4(<pa427a)2) 

xh2(I>J*%i2  +  +  4(<pa424)2) 


(28) 


The  potential  loss  due  to  electron  motion  within  the  anode 
diffuser  is  determined  by  Ohm’ s  law  as 

7 

ohmic  loss  (electronic)  of  anode  diffuser  =  —  (29) 

Aa 

The  diffusion  overpotential  of  the  anode  catalyst  layer  re¬ 
flects  the  potential  loss  caused  by  non-uniform  distribution 
of  hydrogen  concentration  inside  the  layer,  which  is  evalu¬ 
ated  as 


diffusion  overpotential  of  anode  catalyst  layer 
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The  applicability  of  Eq.  (14)  is  limited  to  the  condition: 

.b  \  vAim,a 
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As  Eq.  (31)  is  violated,  Eq.  (16)  is  no  longer  applicable  and 
Eq.  (25)  is  not  valid  either.  Under  such  a  condition,  it  is 
proposed  to  employ  a  piecewise  fourth-degree  polynomial 
for  the  dimensionless  hydrogen  concentration  profile  within 
the  anode  catalyst  layer,  which  is 


4.  Performance  equations  of  a  single  cell 

The  output  potential  of  a  single  cell  can  be  estimated 
as  the  equilibrium  potential  deducted  by  the  potential 
losses  of  its  cathode,  anode  and  membrane  electrolyte. 
Based  on  the  occurrence  of  depletion  of  oxygen  and  hy¬ 
drogen  in  the  catalyst  layers  of  the  cathode  and  the  an¬ 
ode,  four  equations  in  different  regimes  are  formulated 
contingent  upon  the  sign  of  g  and  h,  which  are  defined 
as: 
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where  w  is  the  dimensionless  position  in  the  anode  cata¬ 
lyst  layer  where  hydrogen  depletion  takes  place.  Relating 
the  hydrogen  transport  rate  at  the  interface  of  the  anode 
catalyst-layer  and  the  anode  diffuser  to  the  electric  current 
density  leads  to 
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Accordingly 
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By  following  the  same  procedures,  the  anode  potential  can 
be  expressed  as 


Va  = 


Va°  +  I  1  - 


a-;,  /a 


H2 


I  I 
+ 


a 


A 


a 


1 

+  —  In 


/; 


a 


3/2  +  ,/9/4  +  64^/4(x^)3 


8<Pa/a(*H2)3/2 


(36) 


The  diffusion  overpotential  and  the  electronic  ohmic  loss 
of  the  anode  diffuser  and  the  activation  overpotential  for  the 
hydrogen  oxidation  reaction  have  the  same  expressions  as 
the  previous  condition.  Other  contributive  terms  have  the 
forms 


diffusion  overpotential  of  anode  catalyst  layer 
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For  g  >  0  and  h  >  0: 
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For  g  <  0  and  h  >  0: 
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For  g  <  0  and  h  <  0: 


For  g  >  0  and  h  <  0 
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In  the  above  expressions,  the  potential  loss  of  the  cath¬ 
ode  was  adopted  from  the  results  of  the  previous  work  [27] . 
The  notions  were  modified  so  as  to  be  consistent  with  the 
present  work.  It  should  be  stated  that  the  anodic  reaction  at 
the  cathode  is  neglected  in  formulating  the  cathode  perfor¬ 
mance  equations.  Nevertheless,  both  the  anodic  and  cathodic 
reactions  are  taken  into  account  in  deriving  the  anode  perfor¬ 
mance  equations.  Thus,  it  yields  a  more  complex  form  than 
its  counterpart  for  the  cathode.  Such  an  inclusion  is  compul¬ 
sory  to  provide  accurate  estimations  for  anode  performance 
due  to  a  much  higher  value  of  the  anode  exchange  current 
density.  On  the  other  hand,  the  impact  of  neglecting  the  an¬ 
odic  reaction  at  the  cathode  only  appears  in  an  extremely 
narrow  region  near  the  equilibrium  potential. 


5.  Results  and  discussion 

Because  the  accuracy  of  the  performance  equations  for 
cathodes  has  been  investigated  in  the  previous  study,  this 
study  first  explores  the  differences  between  the  results  cal¬ 
culated  based  on  the  anode  performance  equations  and  the 
corresponding  one-dimensional  model  and  then  discusses 
the  features  exhibited  by  the  performance  equations  of  sin¬ 
gle  cells.  Because  the  chemical  and  physical  properties  of 
PEFC  anodes  might  vary  by  orders  of  magnitude  with  the 
materials,  quantities  and  methods  used  for  fabrication,  the 
validity  of  the  performance  equations  has  to  be  tested  over 
a  wide  range  of  parameter  values.  In  order  to  do  so,  the  pa¬ 
rameter  values  listed  in  Table  1  were  used  as  a  base  case  for 
the  anode,  and  computations  were  carried  out  by  adjusting 


Table  1 


Values  of  physical  and  chemical  parameters  for  anode  of  the  base  case 


Effective  ionic  conductivity  of  ionomer  in 

3  x  10“2 

anode,  kff  (S  cm-1) 

Effective  electric  conductivity  of  anode  gas 

0.5 

diffuser,  crff  (Scm-1) 

Effective  porosity  of  anode  gas  diffuser,  sa 

0.4 

Effective  diffusivity  of  dissolved  hydrogen 

8.4  x 

io-7 

in  the  anode  catalyst  layer,  Dj3]7  (cm2  s-1) 

Effective  gas-pair  diffusivity, 

7.08  x  10"2 

Dh“-C02(=  «a'5DH2-C02)  (cm2  S-1) 

Effective  gas-pair  diffusivity, 

2.11  x  10“2 

£co2-w(=  ea‘5Dco2-w)  (cm2  s-1) 

Effective  gas-pair  diffusivity, 

1.05  x  10-1 

0^-w(=  4'5°h2-w)  (cm2  s-1) 

Thickness  of  anode  gas-diffuser,  da  (fxm) 

300 

Thickness  of  anode  catalyst-layer,  8a  (|xm) 

5 

Product  of  platinum  surface  area  and 

1.4  x 

105 

reference  exchange  current  density  of 

anode,  Aai ®  (A cm-3) 

Anode  pressure,  Pa  (atm.) 

3 

Temperature,  T  (K) 

353 

Bulk  hydrogen  mole  fraction, 

0.633 

Bulk  carbon  dioxide  mole  fraction,  xbCOi 

0.211 

Mole  fraction  of  water  vapor  in  anode,  xWja 

0.156 

Henry’s  constant  for  hydrogen,  H h0 

4.5  x 

104 

(atm.  cm3  per  mole) 

Equilibrium  potential  of  anode,  V®  (V) 

0 

Anode  transfer  coefficient,  aa 

0.5 

the  magnitude  of  a  particular  parameter  with  other  parame¬ 
ter  values  remaining  fixed. 

The  polarization  curves  calculated  using  the  performance 
equations  and  the  one-dimensional  model  for  the  base  case 
are  illustrated  in  Fig.  1.  These  two  curves  almost  coin¬ 
cide.  For  the  parameter  values  used  in  the  calculations,  the 
limiting  current  occurs  around  54  A  cm-2,  which  is  not  in¬ 
cluded  in  Fig.  1  since  such  a  value  is  far  beyond  the  ones 
in  practical  operation.  It  appears  that  the  anode  potential 


Fig.  1.  Anode  polarization  curves  for  base  case:  (•)  one-dimensional 
model;  (■)  performance  equations. 
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Fig.  2.  Polarization  curves  of  individual  contributive  terms  for  base  case 
of  anode  calculated  using  performance  equations:  (■)  activation  overpo¬ 
tential;  (A)  diffusion  overpotential  of  catalyst  layer;  (•)  diffusion  over¬ 
potential  of  gas  diffuser;  (♦)  ohmic  potential  loss  of  ionomer  in  catalyst 
layer;  (T)  ohmic  potential  loss  of  gas  diffuser. 

responds  to  current  density  in  an  approximately  linear  fash¬ 
ion  for  the  calculated  current  range.  Detailed  analyses  of 
the  individual  contributive  terms  of  anode  polarization  are 
displayed  in  Fig.  2.  It  is  seen  that  activation  overpotential  is 
almost  negligible  for  the  current  range  calculated.  This  fact 
is  attributed  to  extremely  fast  reaction  rates  for  hydrogen 
oxidation,  exhibited  by  the  product  value  of  Aa/^  listed  in 
Table  1.  The  ohmic  loss  in  the  anode  diffuser  is  linear  with 
current  density  as  revealed  by  Ohm’ s  law.  It  is  also  seen  that 
the  ohmic  loss  caused  by  proton  migration  within  the  cat¬ 
alyst  layer  increases  in  a  nearly  linear  manner  with  current 
density,  although  its  expression  in  the  performance  equa¬ 
tions  involves  complex  forms  of  fractional  and  exponential 
functions.  For  the  current  range  investigated  in  the  present 
case,  the  diffusion  overpotential  of  the  anode  diffuser  is 
negligibly  small  due  to  relatively  high  values  of  gas-pair 
diffusivity  of  hydrogen  and  the  effective  porosity  of  the 
diffuser.  Compared  with  other  potential  losses,  the  diffusion 
overpotential  of  the  anode  catalyst  layer  varies  with  cur¬ 
rent  density  in  a  more  complex  manner,  which  depends  on 
the  evolution  of  the  hydrogen  concentration  profile  within 
the  layer. 

Various  techniques  have  been  developed  for  the  fab¬ 
rication  of  catalyst  layers  so  as  to  achieve  an  optimum 
structure,  and  hence  improved  performance.  Owing  to  the 
different  nature  of  fabrication  methods,  the  thickness  of  an 
anode  catalyst  layer  may  vary  by  orders  of  magnitude.  As 
a  method  of  deposited  sputtering  is  employed,  a  thickness 
of  less  than  1  [xm  is  achievable  [31].  The  layer  formed  is 
much  thicker,  however,  if  it  is  prepared  by  printing  with 
an  ionomer/carbon-supported-platinum  ink.  In  order  to  take 
such  great  discrepancies  into  account,  two  thicknesses  with 
a  wide  variation  were  selected  for  the  anode  catalyst  layer. 
The  results  are  displayed  in  Fig.  3.  As  shown,  the  polariza¬ 


Fig.  3.  Polarization  curves  for  anodes  with  different  catalyst-layer  thick¬ 
ness.  (A)  8a  =  20  |xm;  (B)  8a  =  0.5  |xm;  (•)  one-dimensional  model; 
(■)  performance  equations. 

tion  curves  obtained  using  the  performance  equations  agree 
well  with  those  from  the  one- dimensional  model.  A  thicker 
catalyst  layer  implies  that  a  larger  space  is  available  to  ac¬ 
commodate  the  hydrogen  oxidation  reaction;  thus,  a  lower 
activation  overpotential  is  expected.  For  an  extremely  fast 
reaction  like  the  present  one,  however,  such  an  effect  is  not 
pronounced.  On  the  other  hand,  a  thicker  catalyst  layer  im¬ 
poses  greater  limitations  for  proton  migration  and  hydrogen 
transport,  which  leads  to  a  greater  loss  of  anode  potential. 
As  illustrated,  the  anode  with  <5a  =  0.5  pirn  exhibits  better 
performance  for  one  with  5a  =  20  pm  for  entire  range  of 
current  density  calculated. 

Investigation  of  the  effects  of  Aa/^  on  the  accuracy  of 
the  present  anode  performance  equations  was  carried  out 
for  two  anodes  with  different  product  values  of  i.e. 
1.4  x  106  and  1.4  x  104  A  cm  3.  The  results  are  presented 
in  Fig.  4.  A  rather  good  agreement  is  also  found  between 
the  performance  equations  and  the  one-dimensional  model, 
as  in  previous  cases.  In  addition,  a  significant  increase  in 
the  anode  potential  loss  is  observed  for  an  electrode  with 
Aa/°  =  1.4  x  104  A  cm  3.  In  order  to  explore  further  the 
sources  of  such  degradation,  individual  potential  losses  rel¬ 
evant  to  Aa/^  were  calculated  based  on  the  performance 
equations.  The  results  are  plotted  in  Fig.  5  and  include  the 
activation  overpotential  for  the  hydrogen  oxidation  reaction 
as  well  as  the  diffusion  overpotentials  of  the  diffuser  and 
of  the  catalyst  layer.  It  is  seen  that  the  diffusion  overpoten¬ 
tials  of  the  gas  diffuser  for  these  two  cases  are  negligible 
and  the  activation  overpotentials  are  only  slightly  different. 
Major  degradation  comes  from  the  diffusion  overpotential 
of  the  catalyst  layer.  Such  results  are  quite  different  to  those 
for  cathodes,  in  which  Acz[?  only  appears  in  the  expression 
of  activation  overpotential  and  the  diffusion  overpotential 
of  the  catalyst  layer  does  not  rely  on  it.  This  dissimilarity 
stems  from  much  higher  values  of  anode  exchange  current 
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Fig.  4.  Polarization  curves  for  anodes  with  different  product  values  of 
Aai°a  (A)  Aai°a  =  1.4  x  104  A  cm-3  (B)  Aai°a  =  1.4  x  106  A  cm-3;  (•) 
one-dimensional  model;  (■)  performance  equations. 

density,  which,  in  turn,  lead  to  a  substantial  participation  of 
the  cathodic  reaction  at  the  anode. 

Among  all  physical  parameters  involved  in  the  perfor¬ 
mance  equations,  is  the  only  one  that  provides  a  quanti¬ 
tative  measure  for  the  transport  ability  of  hydrogen  under  the 
multiphase  configuration  of  the  anode  catalyst  layer.  Since 
the  ionomers  mixed  in  the  catalyst  layer  are  composed  of  hy¬ 
drophobic  backbone  and  hydrophilic  groups,  small  amounts 
of  gas  pores  might  exist,  which  greatly  facilitate  hydrogen 
transport  within  the  layer.  Thus,  the  value  of  may  vary 
greatly  from  one  anode  to  another  due  to  different  degrees  of 
involvement  of  gaseous  hydrogen  transport  within  the  anode 
catalyst  layer.  In  order  to  cover  such  wide  variations,  two 
values  differing  by  a  factor  of  100,  namely  4.2  x  10-7 


Fig.  5.  Polarization  curves  of  individual  contributive  terms  calculated  using 
performance  equations  for  case  of  Fig.  4.  (A)  Aaia  =  1.4  x  104Acm  3; 
(B)  Aai°  =  1.4  x  106  A  cm  3;  (•)  activation  overpotential;  (A)  diffusion 
overpotential  of  catalyst  layer;  (■)  diffusion  overpotential  of  gas  diffuser. 


Fig.  6.  Polarization  curves  for  anodes  with  different  values  of  effective 
hydrogen  diffusivity  in  catalyst  layer.  (A)  =  4.2  x  10  7  cm2  s  1 ;  (B) 

Djk  =  4.2  x  10~5  cm2  s-1;  (•)  one-dimensional  model;  (■)  performance 
equations. 

and  4.2  x  10~5  cm2  s-1  were  employed  for  calculations  and 
the  results  are  given  in  Fig.  6.  As  seen,  quite  good  agree¬ 
ment  is  also  found  between  the  results  from  the  performance 
equations  and  from  the  one-dimensional  model.  A  higher 
suggests  that  hydrogen  is  able  to  diffuse  into  the  in¬ 
ner  part  of  the  catalyst  layer  and  thus  gives  rise  to  increased 
catalyst  utilization.  Consequently,  a  lower  potential  is  ob¬ 
served.  Therefore,  for  the  same  current  density,  the  anode 
with  D^2  =  4.2  x  10  5  cm2  s  1  always  outperforms  that 

with  D^2  =  4.2  x  10-7  cm2  s-1. 

Polarization  curves  for  two  anodes  with  different  values 
of  effective  ionomer  conductivity  in  the  anode  catalyst  layer 
are  shown  in  Fig.  7.  As  illustrated,  the  polarization  curves 


Fig.  7.  Polarization  curves  for  anodes  with  different  values  of  ef¬ 
fective  proton  conductivity  of  ionomer  phase  in  catalyst  layer.  (A) 
ka  l  =  0.01  S  cm-1;  (B)  kff  =  0.1  S  cm-1;  (•)  one-dimensional  model; 
(■)  performance  equations. 
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Fig.  8.  Polarization  curves  for  anodes  with  different  effective  porosity  of 
gas  diffuser.  (A)  ea  =  0.08;  (B)  sa  =  0.5;  (•)  one-dimensional  model; 
(■)  performance  equations. 

based  on  these  two  models  nearly  overlap.  Moreover,  the 
anode  with  kff  =  0.1  Scm-1  always  demonstrates  a  better 
performance  than  that  with  =  0.01  S  cm-1  for  the  same 
current  density.  This  effect  can  be  quite  clearly  revealed 
by  Eqs.  (27)  and  (38)  since  kff  only  appears  in  the  term 
involving  the  ohmic  loss  caused  by  proton  migration  in  the 
performance  equations. 

Due  to  relatively  large  values  of  gas-pair  diffusivity  for 
hydrogen,  the  influence  of  the  effective  porosity  of  the  dif¬ 
fuser  on  anode  performance  only  appears  in  the  range  of  high 
current  densities,  which  are  far  beyond  the  values  of  normal 
operation.  Since  the  objective  of  the  present  study  is  to  exam¬ 
ine  the  accuracy  of  the  performance  equations  by  consider¬ 
ing  the  numerical  solutions  of  the  original  one-dimensional 
model  as  exact  ones,  it  is  also  interesting  to  investigate  the 
accuracy  of  the  performance  equations  in  the  region  near  the 
limiting  current.  In  order  to  achieve  this  purpose,  a  small 
value  of  £a,  namely  0.08,  was  selected  for  calculations  so  as 
to  lower  the  limiting  value  of  current  density.  In  addition, 
another  £a  value  larger  than  that  in  the  base  case  was  also 
chosen  for  comparisons.  The  results  of  both  cases  are  shown 
in  Fig.  8.  As  seen,  the  present  two  models  yield  almost  same 
values  of  anode  polarization.  At  low  current  densities,  where 
diffusion  limitations  of  hydrogen  within  the  diffuser  are  not 
significant,  the  performances  of  these  two  anodes  are  quite 
close.  Nevertheless,  as  the  current  density  is  increased,  hy¬ 
drogen  transport  resistances  within  the  diffuser  play  more 
important  roles  in  anode  performance,  and  thus  the  anode 
with  £a  =  0.08  yields  a  higher  potential  loss  than  the  one 
with  £a  =  0.5.  It  is  also  noted  that  £a  significantly  influ¬ 
ences  the  value  of  the  limiting  current,  which  can  be  exactly 
calculated  by  using  Eq.  (18). 

Comparisons  between  the  polarization  curves  obtained 
from  the  performance  equations  of  a  complete  PEFC  and 
from  the  corresponding  one-dimensional  model  have  been 
carried  out  and  the  results  are  shown  in  Fig.  9.  The  po¬ 


Fig.  9.  Potential  losses  of  cathode,  anode,  membrane  electrolyte  and  single 
cell  for  base  case  (Tables  1  and  2).  (•)  One-dimensional  model;  (■) 
performance  equations. 

tential  losses  pertinent  to  its  cathode,  anode  and  membrane 
electrolyte,  calculated  using  the  performance  equations,  are 
also  illustrated  for  detailed  analysis.  The  values  concern¬ 
ing  the  physical  and  chemical  properties  of  the  cathode  and 
membrane  electrolyte  used  for  the  computations  are  listed 
in  Table  2.  The  effective  porosity  of  the  cathode  diffuser 

Table  2 


Values  of  physical  and  chemical  parameters  for  cathode  and  membrane 
electrolyte  used  for  calculations 


Effective  ionic  conductivity  of  ionomer  in 

3  x  10- 

-2 

cathode,  kff  (Scm-1) 

Effective  electric  conductivity  of  cathode 

0.5 

gas  diffuser,  ofn  (Scm-1) 

Effective  porosity  of  cathode  gas  diffuser,  sc 

0.12 

Effective  diffusivity  of  dissolved  oxygen  in 

5  x  10‘ 

-7 

the  cathode  catalyst  layer,  Dq(  (cm2  s-1) 

Effective  gas-pair  diffusivity, 

2.32  x 

10“3 

£>02-n2(=  £c'5£>o2-n2)  (cm2  s_1) 

Effective  gas-pair  diffusivity, 

3.22  x 

10“3 

£>N2— w(=  £c5d n2-w)  (cm2  s-1) 

Effective  gas-pair  diffusivity, 

3.08  x 

10“3 

£>o2-w(=  £c‘5£o2-w)  (cm2  s-1) 

Thickness  of  cathode  gas-diffuser,  dc  (pan) 

300 

Thickness  of  cathode  catalyst-layer,  8C  (|am) 

10 

Product  of  platinum  surface  area  and 

5  x  10- 

-4 

reference  exchange  current  density  of 

cathode,  Aci°  (A cm-3) 

Cathode  pressure,  Pc  (atm.) 

5 

Bulk  oxygen  mole  fraction,  Vqo 

Bulk  nitrogen  mole  fraction, 

0.190 

0.716 

Mole  fraction  of  water  vapor  in  cathode, 

0.094 

•^w,c 

Henry’s  constant  for  hydrogen,  Ho 2 

2  x  105 

(atm.  cm3  per  mole) 

Equilibrium  potential  of  cathode,  Vc0  (V) 

1.2 

Cathode  transfer  coefficient,  ac 

0.5 

Membrane  thickness,  <5m  (|xm) 

125 

Ionic  conductivity  of  membrane,  km 

0.07 

(S  cm-1) 
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employed  for  calculations  is  smaller  than  its  counterpart  of 
the  anode  so  as  to  account  for  the  effects  of  water  accumu¬ 
lation  within  the  cathode  diffuser.  As  seen,  the  polarization 
curves  constructed  by  these  two  models  are  well  matched. 
According  to  the  results  presented  above  and  in  the  previous 
work,  at  low  current  densities,  a  major  part  of  the  polariza¬ 
tion  comes  from  the  activation  overpotential  for  the  oxygen 
reduction  reaction.  As  the  current  density  is  increased,  over¬ 
potentials  from  other  sources  such  as  ohmic  losses  due  to 
electron  conduction  and  proton  migration,  diffusion  overpo¬ 
tentials  of  the  catalyst  layers  and  the  cathode  gas-diffuser, 
play  more  important  roles  in  cell  performance.  Bear  in  mind 
that  the  activation  overpotential  for  the  hydrogen  oxidation 
reaction  is  always  negligible  due  to  its  extremely  fast  reac¬ 
tion  rates.  In  addition,  owing  to  relatively  large  values  of 
hydrogen  gas-pair  diffusivity,  the  diffusion  overpotential  of 
the  anode  diffuser  is  also  small  for  the  current  range  calcu¬ 
lated.  The  ohmic  loss  of  membrane  is  described  by  Ohm’s 
law,  which  presents  a  linear  dependence  on  current  density. 
The  polarization  curves  shown  in  Fig.  9  reveal  these  conse¬ 
quences.  At  low  current  densities,  cathode  loss  is  predom¬ 
inant.  As  the  current  density  is  increased,  the  anode  and 
membrane  losses  gain  more  importance,  and  substantially 
engage  in  the  overall  cell  overpotential.  Finally,  the  limiting 
value  of  current  density  is  determined  by  the  mass  trans¬ 
port  limitations  of  the  cathode  diffuser  alone.  It  should  be 
pointed  out  that  it  is  not  necessary  to  investigate  the  accuracy 
of  the  performance  equations  for  a  complete  PEFC  since  its 
overall  potential  loss  is  only  a  linear  addition  of  individual 
losses  of  its  anode,  cathode  and  membrane  electrolyte.  The 
accuracy  of  the  performance  equations  for  cathode  and  an¬ 
ode  has  been  investigated  in  the  present  and  previous  stud¬ 
ies,  in  addition,  the  membrane  loss  can  be  estimated  exactly 
because  of  its  linearity.  Therefore,  a  high  accuracy  of  the 
performance  equations  for  a  complete  PEFC  is  concluded 
for  the  parameter  range  investigated. 

6.  Summary  and  conclusions 

A  one-dimensional  model  has  been  formulated  to  take 
into  account  the  transport  and  electrochemical  processes 
within  the  anode  catalyst  layer  and  diffuser  of  PEFCs, 
and  to  produce  performance  equations  in  algebraic  form. 
The  reduction  procedures  are  performed  by  employing  a 
fourth-degree  polynomial  approximation  or  a  piecewise  one 
for  the  hydrogen  concentration  profile  within  the  catalyst 
layer  contingent  upon  the  occurrence  of  hydrogen  deple¬ 
tion  at  the  interface  between  the  membrane  and  the  anode 
catalyst-layer.  After  lumping  the  reaction  rate  at  a  reaction 
centre,  the  performance  equations  in  algebraic  form  are  de¬ 
rived.  Like  the  equations  for  cathodes  developed  previously, 
individual  overpotential  terms,  resulting  from  the  limiting 
processes  considered  in  the  original  one-dimensional  model, 
can  be  quantitatively  estimated  for  the  anode.  Computa¬ 
tional  results  show  that  the  polarization  curves  created  by  the 


anode  performance  equations  and  by  the  one-dimensional 
model  agree  well  for  the  extensive  parameter  range  inves¬ 
tigated.  The  performance  equations  for  a  complete  PEFC 
are  formulated  by  linearly  combining  the  potential  losses 
pertinent  to  the  cathode,  anode  and  membrane  electrolyte. 
The  characteristics  of  the  polarization  curves  of  single  cells 
have  been  discussed.  It  is  shown  that  although  the  anode 
activation  overpotential  can  be  always  neglected  due  to  the 
fast  rates  of  the  hydrogen  oxidation  reaction,  other  anode 
potential  losses,  such  as  diffusion  overpotential  of  the  cata¬ 
lyst  layer  and  ohmic  losses  caused  by  proton  migration  and 
electron  conduction  are  not  negligible  at  medium  and  high 
current  densities.  Thus,  the  present  performance  equations 
that  include  the  potential  losses  of  the  cathode,  anode  and 
membrane  electrolyte  are  able  to  provide  accurate  evalua¬ 
tion  of  the  performance  of  a  complete  PEFC. 
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